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It is well known that the [O iii]/l5007 emission line and 
hard X-ray(2-10keV) luminosities are good indicators of AGN 
activities and that the near and mid-infrared emission of AGN 
originates from re-radiation of dusty clouds heated by the 
UV/optical radiation from the accretion disk. In this paper 
we present a study of the near-infrared K-band (2.2yum) prop- 
erties for a sample of 65 Seyfert 2 galaxies. By using the 
AGN/Bulge/Disk decomposition technique, we analyzed the 
2MASS Ks-band images for Seyfert 2 galaxies in order to de- 
rive the Ks-band magnitudes for the central engine, bulge, and 
disk components. We find that the Ks-band magnitudes of the 
central AGN component in Seyfert 2 galaxies are tightly cor- 
. . . , related with the [O iii]/15007 and the hard X-ray luminosities, 
■ which suggests that the AGN K-band emission is also an ex- 
, ] cellent indicator of the nuclear activities at least for Seyfert 2 
■ " galaxies. We also confirm the good relation between the central 
, . \ black hole masses and bulge's K-band magnitudes for Seyfert 
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1. Introduction 

In the standard unified scheme, Seyfert 1 and 2 galaxies are in- 
trinsically the same objects; the absence of broad Balmer emis- 
sion lines in Seyfert 2s is due to the obscuration by a pc-scale 
dusty torus oriented along the line of sight (Krolik & Begelman 
1988; Antonucci 1993). The primary evidence for this model 
is the detection of polarized broad Balmer emission lines in 
dozens of Seyfert 2 galaxies, which have betrayed the hidden 
broad line regions via scattering off electrons and/or dust lo- 
cated above the hole of the torus (Antonucci & Miller 1985; 
Tran 1995 and 2001; Young et al. 1996; Heisler et al. 1997; 
and Moran et al. 2000). 

One natural prediction of the unified model is that such 
a dusty torus surrounding the central engine will absorb the 
UV/optical radiation from the accretion disk and re-radiate 
it in the infrared, as the sublimation temperature of graphite 
grains is in the range 1500 - 1800 K. According to the Wien 
law, TApeakij^Jn) = 2898, the peak wavelength of re-radiation 
from the dusty torus should be around 2//m. In fact, it is well 
known that Seyfert galaxies are strong near and mid-infrared 
sources (Rieke 1978; Neugebauer et al. 1979). Pier & Krolik 
(1992) performed the first calculation of thermally re-radiated 
infrared spectra of the compact dust torus, and found good 
agreement between the model predictions and observations in 
Pier & Krolik (1993), which suggested that near infrared emis- 
sion in Seyfert galaxies mainly arise from the re-radiation of 
dusty torus heated by UV/optical light from the central engine. 
Recently, Suganuma et al. (2005) presented reverberation mea- 
surements of several well-known Seyfert 1 galaxies by mon- 
itoring flux variations in the optical (UBV) and near-infrared 
(JHK) bands. They find a clear time-delayed behavior of the K- 
band flux variations relative to the V-band. From the flux vari- 
ation gradients they derive (H - K) color temperatures of 1500 



Send offprint requests to: Q. Gu 

* Visiting Scholar, Harvard-Smithsonian Center for Astrophysics, 
60 Garden St., Cambridge, MA 02138, USA 



- 1800 K, and suggest that the bulk of the K-band emission 
originates from thermal emission of dust grains in the optically 
thick torus (see also Glass 2004). Thus the K-band emission 
should also be a good indicator of AGN activity. 

The recent, comprehensive Two Micron All Sky Survey 
(2MASS) has uniformly scanned the entire sky in three near- 
infrared (NIR) (J 1.25//m, H 1.65//m, and Ks 2.11 pm) bands. 
The point-source sensitivity limits (lOcr) are J= 15.8 (0.8 mJy), 
H=15.1 (1.0 mJy), and Ks = U3 (1.3 mJy) mag. The ex- 
tended source sensitivity (lOcr) is ~ 1 mag brighter than the 
point source limits^ The 2MASS public catalog contains > 
1,000,000 galaxies and thus provides the most comprehensive 
database for the study of NIR galaxy properties. 

In this paper we have performed AGN/Bulge/Disk decom- 
position for a sample of Seyfert 2 galaxies by using 2MASS 
Ks-band images in order to derive the near-infrared magnitudes 
of the central AGNs. In Sect. 2, we present a short descrip- 
tion for our sample and present the results of decomposition 
in Sect. 3. We discuss our results in Sect. 4 and present conclu- 
sions in Sect. 5. By using the WMAP data, Spergel et al. (2003) 
have well determined the Hubble constant to be Ho= 71 km s~^ 
Mpc~\ which will be used throughout this paper. 

2. Data 

We recently performed a systematic study of a large and ho- 
mogeneous sample of 65 nearby Seyfert 2 galaxies observ- 
able from the southern hemisphere. A full description of the 
sample selection, observation, and data reduction is presented 
in Joguet et al. (2001). The star formation history and stellar 
populations in the central ~ 200pc for these Seyfert 2s are de- 
rived by means of stellar population synthesis modelling (Cid 
Fernandes et al. 2004), and Gu et al. (2005) analyzed the "pure" 
emission line spectra obtained by subtracting the synthetic stel- 
lar contribution to derive pure emission-line spectra. Here we 
use the 2MASS Ks-band images of all Seyfert 2 galaxies in our 
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sample^, to present a study of NIR properties for these galaxies. 
The plate scale of the 2MASS Ks-band image is 1.0 arc-second 
per pixel. 

2.1. Light decompositions 

We applied the 2-dimensional bulge-disk decomposition pro- 
gram, GALFIT VERSION 2.0.3c (Peng et al. 2002), to the 
2MASS Ks-band images of these 65 Seyfert 2 galaxies. 
Specifically, in order to separate the AGN component, we fit- 
ted each galaxy with a combination of a point spread func- 
tion (PSF) centered in the image, an inclined exponential disk, 
and a bulge with a surface brightness profile of ^x;?[-(r/r^)^/"] 
(Sersic 1968) with the index n constrained to the range of 0.5 
<n<5.0. 

Ground-based images usually show near-Gaussian or 
Moff'at-like PSFs (Peng et al. 2002). We generated an analytic 
profile of PSF for the 2MASS Ks-band images using the stan- 
dard tasks in IRAF^. Following Jarrett et al. (2000), we used 
a generalized, radially symmetric exponential function to de- 
scribe the 2MASS Ks-band PSF, which is: 



/(r) = /oexp 



(1) 



where /o is the central surface brightness, and a and JS are free 
parameters. For the case of the 2MASS survey, we adopted JS = 
0.5 and a = ^f2o- = 1.5'^ where cr = FWHM/2.354, FWHM 
is the seeing of the 2MASS images, and FWHM= 2.5'' for the 
Ks images ^. 

Examples of the fitting results are shown in Fig.l for 4 rep- 
resentative Seyfert 2s covering the full range of;^^ (NGC 3035, 
NGC 5252, NGC 5643, and NGC 7582 with of 0.012, 0.013, 
0.032, and 0.049, respectively). The standard measurement er- 
rors for AGN, bulge, and disk Ks-band magnitudes are 0.2, 
0.3, and 0.5 mag, respectively. For NGC 5643 there appears 
to be a polar ring left after subtracting the bulge/disk compo- 
nents. This illustrates the high quality of the Bulge/Disk sub- 
traction and may indicate that the ring reflects a dynamical in- 
teraction with another galaxy. The currently favored model for 
the formation of polar rings is suggested to be galactic mergers 
(Sparke & Cox 2000; Bournaud & Combes 2003; and Maccio, 
Moore & Stadel 2006). For NGC 7582, the residual flux after 
B/D subtraction is peculiar, which may indicate a warp or some 
other faint asymmetry in the disk. It is interesting to note that 
for three objects (Mrk 897, Mrkl370, and NGC 788 with mor- 
phological types of Scd, Sa, and SA(s)0/a, respectively), there 
seems to be no need for exponential disk components to fit the 
Ks-band images. 

The decomposed Ks-band magnitudes of AGN, bulge, and 
disk are listed in Table 1 , where for completeness we also tab- 
ulate the stellar velocity dispersions (cr) from Cid Fernandes et 
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Fig. 2. The AGN Ks-band absolute magnitude versus [O 
iii]/15007 luminosity. Red open-square symbols correspond to 
Seyfert 2s with direct evidence of hidden broad-line regions 
(HBLR) from spectropolarimetric observations. The dashed 
line shows a ordinary least square (OLS) bisector fit to the data 
of slope -0.715 ± 0.051. The Spearman rank-order correlation 
coefficient (Rs) is -0.52 and the nuU probability is 0.0004. 

al. (2004), and the [O iii]/l5007 emission line and the hard X- 
ray(2-10 keV) continuum luminosities from Gu et al. (2005). 

3. Results 

In the past three decades many attempts were made to es- 
tablish which are the truly isotropic emissions from Seyfert 
galaxies. Among various emissions at diff'erent wavelengths, 
the [O iii]/l5007 and hard X-ray (2-lOkev) continuum luminosi- 
ties have been found to show similar distributions for Seyfert 1 
and 2 galaxies (Dahari & De Robertis 1988; Keel et al. 1994; 
Mulchaey et al. 1994; Alonso-Herrero et al. 1997; Kauff'mann 
et al. 2003), indicating isotropic emission and thus making 
them reliable tracers of the intrinsic nuclear activity. Alonso- 
Herrero, Ward & Kotilainen (1997) compared the hard X-ray, 
the [O iii]/15007, and the small-aperture L-band (3.5yum) prop- 
erties of low-redshift PG quasars, Seyfert 1, and 2 galaxies, and 
they found good correlations between the [O iii]/l5007 and L- 
band properties and between the hard X-ray and L-band prop- 
erties for Seyfert Is and PG quasars. 

Figure 2 shows the correlation between the Ks-band abso- 
lute magnitudes of our Seyfert 2s with [O iii]/15007 luminosity. 
The red open squares represent objects where the spectropolari- 
metric observations detected the polarized broad emission lines 
(see also Gu et al. 2005), mainly from the compilation by Gu & 
Huang (2002) plus new observations by Lumsden, Alexander 
& Hough (2004). As both the AGN Ks-band absolute magni- 
tude (Mk,agn) and [O iii]/15007 luminosity (L[oiii]) are well de- 
termined, in order to get the best fit between Mk,agn and L[oiii] , 
we used both the ordinary least- square (OLS) bisector method 
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Fig. 1. Results of decomposition. As an example, we only show 4 representative Seyfert 2s with the smallest to the largest 
(NGC 3035, NGC 5252, NGC 5643 and NGC 7582,from top to bottom). For each source, we show the original image (left), the 
model (middle) fitted by GALFIT, and the residual image (right). 



(Isobe et al. 1990) and the Buckley- James regression method 
from the Astronomy Survival Analysis package (ASURV vl.2; 
Isobe et al. 1985; 1986), and derived the same fitting results, 
which is logL[oiii] = (-0.715 ± 0.051) Mk,agn + (25.260 ± 
1.092). Figure 2 also quotes the Spearman rank-order correla- 
tion coefficient (Rs, Press et al. 1992) as -0.52, and a probabil- 
ity of Pnuii = 0.0004 for the null hypothesis of no correlation 
between Mk,agn and L[oiii]. 



Similarly, Fig. 3 shows the correlation between the AGN 
Ks-band absolute magnitude versus the hard X-ray 2-lOkeV 
luminosity. Twelve of our Seyferts are Compton-thick sources 
with column densities larger than lO^'^cm"^, which are de- 
rived by fitting a simple power law with cold absorption to 
the X-ray spectrum. For Compton-thick objects, the direct hard 
X-ray continua are completely absorbed, so we cannot es- 
timate their intrinsic hard X-ray luminosities. For these ob- 
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Fig. 3. The AGN Ks-band absolute magnitude versus 2-10 keV 
hard X-ray luminosity. The dashed line shows an OLS bisector 
fit to the data of slope -0.406+0. 124. The Spearman rank-order 
correlation coefficient (Rs) is -0.54 and the null probability is 
0.0008. Symbols have the same meaning as in Fig. 2. 

jects, the hard X-ray 2-lOkeV luminosities will be taken as 
lower limits. We have to derive the linear regression using 
the Buckley- James regression method from ASURV, which is 
logL2-iokev = (-0.406+0. 124) Mk,agn+(33.762±2.594), with 
a Spearman rank-order correlation coefficient (Rs) of -0.54 
and a null probability of Pnuii = 0.0008. As a worst case, if 
ALL lower limits were detected at their exact lower limit val- 
ues, the correlation still holds but becomes much steeper. By 
using the same survival techniques, we derive the best fit to be 
logL2-iokeV = (-0.912+0.056) Mk,agn+(22.647+1.237), with 
the Spearman rank-order correlation coefficient (Rs) of -0.30 
and a null probability of Pnuii = 0.061. 

Both Fig. 2 and Fig. 3 indicate that the AGN K-band mag- 
nitudes are tightly correlated with [O iii]/l5007 and the hard 
X-ray(2-10keV) luminosities, both of which are considered to 
be good indicators of nuclear activities. Thus it suggests that 
the AGN K-band magnitude is also a good indicator of AGN 
activities for Seyfert 2 galaxies. We find in Fig. 2 that Seyfert 
2s with hidden broad-line regions clearly have much higher 
[O iii]/l5007 luminosities and higher Ks-band absolute mag- 
nitudes, indicating much powerful AGN activities as found by 
Gu & Huang (2002). 

4. Discussion 

Since the seeing of 2MASS Ks-band images is rather poor 
(~ 2.5''), it is possible that part of the emission from AGN 
is mixed with that of the bulge. In order to check whether 
GALFIT did work in this situation, we used the Hubble Space 
Telescope NICMOS K-band images, which have much better 
spatial resolution than 2MASS. Fortunately, 3 sources in our 
sample have NICMOS K-band data: NGC 2110, NGC 2992, 



and NGC 5643. The AGN magnitudes from NICMOS GALFIT 
fitting are 12.17, 12.03, and 10.96 mag, respectively. Our val- 
ues from the 2MASS data are 12.15, 11.33, and 10.69 mag, re- 
spectively. For NGC 2110 and NGC 5643, the fitting results are 
consistent within the typical fitting errors of GALFIT (see Sect. 
2) with the NICMOS magnitudes. But for NGC 2992, which 
is an unabsorbed Seyfert 2 galaxy with a polarized broad Ha 
emission (Panessa & Bassani 2002; Rix et al. 1990), the magni- 
tude diff'erence is as large as 0.7 mag, which could be evidence 
of variability or perhaps of an aperture/PSF eff'ect. 

4.1. Relation between BH mass and bulge luminosity 

It is well known that there is a tight correlation between 
black hole masses and the bulge properties in spiral galax- 
ies (Magorrian et al. 1998; Gebhardt et al. 2000; Ferrarese & 
Merritt 2000). Marconi & Hunt (2003) find that NIR bulge 
luminosities and BH masses are tightly correlated. For those 
galaxies with accurate BH masses, the spread of such a rela- 
tion is nearly the same as the famous Mbh - o- relation. 

Figure 4 shows the relation of bulge K-band absolute mag- 
nitudes and BH masses for our sample of Seyfert 2s. The BH 
mass is derived from the well-known Mbh - cr relation, which is 
log(MBH/Mo) = (8.13+0.06)+(4.02+0.32)xlog(cr/200kms-^) 
(Tremaine et al. 2002), where cr is the stellar velocity disper- 
sion derived by Cid Fernandes et al. (2004). The dashed line is 
the best fit to the data by means of the OLS bisector method, 

which is logMBH(Mo) = (1.199 + 0.080) log(LK,Bulge/LK,0) - 

(5.184+0.831). It is very interesting to note that the relation be- 
tween bulge K-band magnitude and BH mass for Seyfert 2s is 
nearly the same as that derived for normal galaxies with secure 
BH mass measurement (see Table 2 in Marconi & Hunt 2003), 
which shows that host galaxies of Seyfert 2s obey the same re- 
lation as normal galaxies. However, our scatter is much larger 
than that of Marconi & Hunt (2003), the main reason being 
that BH masses are well-measured in Marconi & Hunt (2003), 
while we estimate BH masses via the Mbh - cr relation, and the 
stellar velocity dispersion was derived by stellar synthesis and 
thus has larger uncertainties. 

5. Conclusion 

In this paper, we present a study of NIR properties for a sample 
of Seyfert 2 galaxies. After decomposing the 2MASS Ks-band 
image in the AGN, bulge, and an exponential disk, we derived 
the Ks-band magnitudes for each component. We find that the 
AGN Ks-band magnitudes are tightly correlated with the lu- 
minosities of [O iii]/l5007 and the hard X-ray emission, which 
suggests that the K-band emission is also an excellent indicator 
of nuclear activity, at least for Seyfert 2 galaxies. We also con- 
firm the good relation between the central black hole masses 
and bulge's K-band magnitudes. 
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9.77 





40. 


,16 




NGC 1068 


16.1 


143.5 


8.02 


7.32 


6.80 


42, 


.00 


>41.03'' 


NGC 1125 


46.2 


104.7 


12.05 


11.38 


10.77 


39. 


,84 




NGC 1667 


64.0 


148.9 


12.20 


9.63 


10.42 


jy. 




>42.39'' 


NGC 1672 


19.0 


96.8 


11.54 


8.95 


7.82 


38. 


,88 


41.21 


NGC 2110 


32.9 


241.5 


12.15 


8.53 


10.02 


39. 


.99 


42.62 


NGC 2979 


38.3 


112.2 


12.62 


10.39 


10.12 


39. 


,29 




NGC 2992 


32.5 


171.8 


11.33 


8.96 


8.47 


40. 


.46 


41.76 


NGC 3035 


61.4 


161.4 


12.06 


12.41 


10.03 


40. 


,41 




NGC 3081 


33.6 


133.4 


12.28 


10.79 


9.63 


40. 


,55 


41.96 
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Table 1. continued 



Zhixin Peng et al.: The K-band properties of Seyfert 2 galaxies 



Galaxy 


D 


cr 


itik,agn 


lTlK,bulge 


mK,disk 


log L[oiii] 


log L2_i0keV 




(Mpc) 


(km s-i) 


(mag) 


(mag) 


(mag) 


(ergs s-^) 


(ergs s~^) 


NGC 3281 


45.1 


160.0 


13.24 


10.22 


8.54 


40.26 


42.84 


NGC 3362 


116.7 


103.5 


13.68 


11.05 


13.19 


41.59 




NGC 3393 


52.8 


156.7 


11.47 


10.60 


9.83 


41.05 


>41.13" 


NGC 3660 


51.8 


95.5 


12.39 


10.51 


11.57 


40.63 




NGC 4388 


35.5 


110.7 


11.26 


9.23 


9.22 


40.98 


42.81 


NGC 4507 


49.9 


144.5 


11.44 


9.51 


10.70 


41.22 


43.33 


NGC 4903 


69.5 


200.4 


13.33 


10.46 


12.41 


41.02 




NGC 4939 


43.8 


154.9 


12.18 


11.50 


9.04 


40.75 


42.01 


NGC 4941 


15.6 


98.2 


11.25 


9.01 


9.22 


39.99 


40.95 


NGC 4968 


41.6 


121.1 


12.05 


10.75 


9.44 


40.49 


>40.92^ 


NGC 5135 


57.9 


142.6 


11.11 


10.17 


9.98 


40.57 


>40.91^ 


NGC 5252 


97.4 


208.9 


12.76 


10.49 


13.44 


40.49 


43.12 


NGC 5427 


36.9 


100.2 


13.16 


9.88 


9.29 


39.90 




NGC 5506 


26.1 


97.9 


9.63 


10.74 


9.02 


40.74 


42.94 


NGC 5643 


16.8 


92.9 


11.17 


7.65 


11.25 


39.91 


> 40.65" 


NGC 5674 


105.2 


128.5 


11.99 


12.27 


10.21 


39.92 


43.23 


NGC 5728 


39.3 


155.2 


12.43 


9.89 


8.41 


40.82 




NGC 5953 


27.7 


92.9 


12.88 


10.03 


10.38 


39.31 




NGC 6221 


20.8 


111.4 


10.72 


8.21 


7.96 


39.33 




NGC 6300 


15.6 


99.8 


11.22 


9.65 


6.15 


39.44 


41.28 


NGC 6890 


34.1 


108.9 


12.40 


9.35 


13.41 


40.64 




NGC 7172 


36.7 


190.0 


12.18 


10.23 


8.89 


39.00 




NGC 7212 


112.5 


167.9 


12.47 


11.70 


9.87 


42.01 


42.29 


NGC 7314 


20.0 


59.9 


12.86 


11.43 


7.71 


39.17 


42.29 


NGC 7496 


23.2 


100.7 


12.15 


9.91 


8.87 


39.38 


41.71 


NGC 7582 


22.2 


132.1 


10.13 


8.98 


7.43 


40.39 


42.21 


NGC 7590 


22.5 


98.6 


11.53 


9.63 


8.44 


39.26 


40.86 


NGC 7679 


72.4 


95.9 


13.63 


10.59 


10.27 


40.31 


42.52 


NGC 7682 


72.4 


152.4 


14.65 


11.21 


9.69 


40.74 


42.91 


NGC 7743 


24.1 


95.3 


11.60 


9.98 


7.97 


38.74 


39.56 



" Compton- thick source with column density > 10^"^ cm ^. 



